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first time which gives a small contribution for light targets but is of the same or- 
der of magnitude as nuclear breakup for heavy targets. The overall interference 
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Thus it appears that Coulomb breakup experiments need to be analyzed with 
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1 Introduction 



This paper is concerned with the nuclear and Coulomb breakup of weakly bound 
nuclei, in particular halo nuclei in which the last neutron is very weakly bound 
and in a state of low angular momentum. As a consequence its distance from 
the core center is quite large and this makes the nuclear and Coulomb breakup 
strong effects in particular when the target is heavy. 

There have been a large number of approaches pj|-|6^| attempting to give 
an accurate description of this phenomenon and the most recent papers contain 
a review of the present situation [j], ^, t|. The accuracy of the method is an 
important issue not only from the point of view of the understanding of the 
reaction mechanism but also because a number of parameters of the neutron 
initial state wave function such as the binding energy and the spectroscopic 
factor can be deduced from the data analysis 

An important aspect of the problem is that nuclear breakup is always present 
at the same time. It is dominant for a light target and it remains non negligible 
for heavy targets. Its contribution to the total breakup cross section on a 
heavy target has been a matter of debate but it should range from about 20% 
to 50%. In such conditions one expects to see interference effects. However 
most of the analysis of experimental data so far performed have relied on the 
assumption that interference effects are small in the intermediate energy range 
(40-80 A.MeV). Furthermore nuclear and Coulomb breakup have been treated 
as independent processes and calculated in the absence of each other to different 
orders in the interactions. Such approximations have been justified by the fact 
that the nuclear interaction is strong and has a short range, therefore needs to be 
treated to all orders. The Coulomb potential on the other hand is comparatively 
weaker and has a long range, therefore it has often been treated to first order. 

One difficulty in the theoretical approach comes from the fact that one needs 
to have a consistent description of the nuclear and Coulomb breakup amplitudes 
from which the interference originates. In this respect the most accurate calcula- 
tions presently available are the numerical solutions of the Schrodinger equation 
H . They treat both interactions to all orders and therefore contain also all pos- 
sible interference effects. However it is not so easy to extract from them a simple 
picture of the coupling mechanisms. The numerical solution of the Schrodinger 
equation has also been used with only the nuclear || or the Coulomb potential 
[ fio[ to test the accuracy of simpler models used to calculate the nuclear 

and Coulomb breakup separately. The conclusions of recent works have been 
that, for well developed haloes, nuclear breakup needs to be treated to all orders 
but the sudden and eikonal approximations are acceptable since they give an 
accuracy of about 20% || [l3| compatible with the present experimental error 
bars. On the other hand it seems that Coulomb breakup can be safely treated 
to first order jy, [|, [l0| , but the time dependence needs to be retained in order 
to reproduce some very fine experimental observables such as the neutron-core 



final energy spectrum. DWBA [LJ, [16|, scmiclassical |17], |18| [19|, adiabatic 
[ pp| and coupled channel (2lJ] approaches have also been proposed. A common 
feature of some of those approaches J2], |l2| 18, [T^] has been to give a smaller 
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Coulomb breakup contribution than that calculated with the first order pertur- 
bation theory. Also many of them were concerned with proton breakup pl| . p2[ 
and lower energy regimes. Such situations are physically more challenging and 
we do not attempt to include them in the present discussion. Therefore we 
restrict ourselves to neutron breakup at intermediate energies in the attempt 
to clarify the role of Coulomb-nuclear couplings and interference and to make a 
link between numerical solutions and simple analytical methods. 

An important issue is to determine the experimental observable on which 
couplings and interference effects would show up best. Interferences are a quan- 
tum mechanical manifestation of a coherence effect. Therefore in order to see 
them an observable of the "exclusive" type is expected to be most suitable. 
Some work has been done jl7|, 23, |24| to study interference effects on the 
ejectile angular distribution but nothing exists at present on the interference 
effects on the neutron angular distribution which we address here. Neutron an- 
gular distributions have been the best experimental proof of the existence of the 
two distinct mechanisms of Coulomb and nuclear diffractive breakup. Further- 
more they give information on the relative strength of the two processes, on the 
target dependence and beam energy dependence of the reaction mechanisms. 

Neutron angular distributions due to nuclear and Coulomb breakup were 
measured and calculated in |25| ] but information on the interferences were not 
extracted. Another reliable calculation of the neutron angular distribution from 
nuclear breakup alone was given in [Q. Models for the reaction mechanism 
and structure information extracted from the data can be considered reliable 
only if three aspects of the data are reproduced at the same time within a 
given model. One is the ejectile or neutron energy spectrum or its equivalent 
parallel momentum distribution. Another is the neutron angular distribution. 
Finally the absolute breakup cross section is important to fix the initial state 
parameters. 

In this paper we will try to get some insight on the problem by a simple but 
consistent calculation of the Coulomb and the nuclear breakup in the presence 
of each other and of their interference. This will allow us to get both the neutron 
angular distribution as well as its energy spectrum due to the combined effects 
of the nuclear and Coulomb interaction. The absolute breakup cross section is 
obtained from the integration of the angular distribution. We will discuss under 
which experimental conditions interference effects could be seen and why they 
have not been clearly identified so far. An important result of our method is 
that it gives simple analytical expressions for the Coulomb cross section for any 
initial angular momentum state of the neutron. 

Section 2 of the paper presents some theoretical background. The purpose is 
to clarify the assumptions made in deriving the nuclear and Coulomb breakup 
amplitudes. Section 3 gives expressions for the cross-section and the breakup 
amplitudes. The results of calculations are presented in section 4 and conclu- 
sions in section 5. Some details of the derivation of the Coulomb amplitude are 
given in the first Appendix. The calculations presented in the paper assume 
an s-state initial neutron wave function in the projectile. Some formulae which 
generalize the results for the Coulomb amplitude are given in the second ap- 
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pcndix. They are included for completeness and are not used in the present 
paper. 

2 Theoretical background 

We consider the breakup of a halo nucleus like 11 Be consisting of a neutron 
bound to a 10 Be core in a collision with a target nucleus. In previous works on 
halo breakup (2(| we have assumed that the projectile core and the target 
move along classical trajectories. Here we adopt a different approach. The 
system of the halo nucleus and the target is described by Jacobi coordinates 
(R, r) where R is the position of the center of mass of the halo nucleus relative 
to the target nucleus and r is the position of the neutron relative to the halo 
core, and the coordinate R is assumed to move on a classical path. This allows 
target recoil to be included in a consistent way. The Hamiltonian of the system 
is 

H = T R + T r + V nc (r) + V 2 (R,r) (1) 

where Tr and T r are the kinetic energy operators associated with the coordinates 
R and r and V nc is a real potential describing the interaction of the neutron 
with the core. The potential V% describes the interaction between the projectile 
and the target. It is a sum of two parts depending on the relative coordinates 
of the neutron and the target and of the core and the target 

V 2 (R, r) = V nt (for + R) + V ct (R - /V) (2) 

Here /3i = m n /m Pl P2 = m c /m p = 1 — j3i, where m n is the neutron mass, m c is 
the mass of the projectile core and m p = m n + m c is the projectile mass. Both 
V n t and V c t are represented by complex optical potentials. The imaginary part 
of V n t describes absorption of the neutron by the target to form a compound 
nucleus. It gives rise to the stripping part of the halo breakup. The imaginary 
part of V c t describes reactions of the halo core with the target. The potential 
V c t also includes the Coulomb interaction between the halo core and the target. 
This part of V c t is responsible for Coulomb breakup. 

The mass ratio Pi is small for a halo nucleus with a heavy core. For example 
(3i sa 0.1 and fo ~ 0.9 in the case of 11 Be. This property is used here to 
approximate the neutron-target and neutron-core potentials by 

V nt (p 2 r + R) « Vnt(r + K) (3) 

K*(R-/V) « K t (R)+V e// (r,R) (4) 

where 

V e //(r,R) = ftr • F rt (R) and F ct (R) = -W ct (R) (5) 
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Here P ci (R) is the classical force acting between the target and the projec- 
tile core. The halo breakup is caused by the direct neutron target interaction 
V n t or by a recoil effect due to the core-target interaction. Coulomb breakup 
of a one-neutron halo nucleus is a recoil effect due the Coulomb component V c t 
of the core-target interaction and is contained in V e // (r, R). ft is proportional 
to the mass ratio (3\. The dipole expansion Eqs.(^),(||) is a good approximation 
provided fi\r << R. In the case considered here r is of the order of the halo 
radius (4 to 6 fm) and R is larger than the core-target strong absorption radius 
( about 6 fm for a 9 Be target and 12 fm for a Pb target). As (3\ » 0.1 for the 
n Be halo nucleus the dipole approximation should be valid. 

The nuclear part of the core-target interaction can also contribute to V e j / 
and gives rise to the 'shake-off' component of the breakup amplitude which is 
discussed in ref. |28|]. They show that it can be important for a light target. It 
is not included in the present paper. 

The theory in this paper is based on a time dependent approach which can 
be derived from an eikonal approximation. The projectile motion relative to 
the target is described by a time-dependent classical trajectory R (t) = d + vtz 
with constant velocity v and impact parameter d ( z is a unit vector parallel 
to the z-axis). With the approximations (||j4|) and (|5|) to the potentials the 
wave function <j> (r, d, t) describing the dynamics of the halo neutron satisfies 
the time-dependent equation 

•<.' Mr ' d ' /: ) = (H r + V nt (r + K(t))+V eff (r,K(t)))ct ) (r,d,t) (6) 



— - yi r -r K„ tV x -r u W , -r v e// 

where H r = T r + V nc (r) is the Hamiltonian for the halo nucleus. As t — + — oo 
the wave function tends to the initial halo nucleus wave-function 

</> (r,d, i) -> (j)i m (r, t) = 4>i m (r) exp (-~ie t/fr) (7) 

with binding energy £o- The initial state (\>i m (r) can be an s-state or a general 
single particle state with angular momentum (l,m). 

In the present paper we neglect the final state interactions between the neu- 
tron and the halo core, but include the final state interactions between the 
neutron and the target. This approximation should be satisfactory unless there 
are resonances in the neutron-core final state interaction which give rise to mod- 
ification in the momentum distributions of the fragments of the type recently 
discussed in J^]. Such effects have not been seen so far in the experimental 
distributions. When the neutron-core final state interactions are neglected the 
breakup amplitude can be written as 

If 00 

g lm (k,d) = - / dt (<f> f (t) \V 2 (r, t) \<f> lm (f)) (8) 



where V2 (r,t) — V nt (r + R (i) ) + V e ^^ (r, R (t)) and (f>f (t) satisfies the equa- 
tion 

in?^P- = (T r +V 2 (r,t))<i> f (t) (9) 
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with the boundary condition that 4>f (t) ~ exp (ik • r — ie^t/K) when t is large. 
The final step is to make an eikonal approximation for (f>j (t) 

4> } (t) = exp(ik- r- ie k t/?i) exp J V 2 (r,t') dt'j (10) 

The breakup amplitude becomes 

g lm (k,d) = 1 /" d 3 r /" dte- 4k - r +^ t e(* T M r .*')*') y 2 ( r , 4 ) <^ m ( r ) (11) 

where w = (e^ — £q) /7i. 

The components V^j and V e // of V2 (r,t) have to be treated differently be- 
cause of the long range of the Coulomb interaction in V e f / . The neutron-target 
interaction is strong and has a short range. In this paper we assume that the 
interaction time tj for this part of the interaction is very short in the sense that 
ojTf is small compared with unity. On the other hand the dominant contribution 
to V e ff comes from the long range Coulomb interaction between the halo core 
and the target. It is weaker and changes more slowly. In this paper we calcu- 
late the contribution of V e ff to first order but include its full time dependence. 
We will show that under such hypothesis the total amplitude gi m (k,d) can be 
written as a sum of two parts 

g lm (k,d) = CT' (k,d) + g? n T d (k,d) (12) 

where g^ cl (k,d) is zero order in V e ff and gf^ ul (k,d) is first order. Explicit 
expressions for the Coulomb and nuclear contributions are given in the next 
section. 



3 Amplitudes and Cross Sections 
3.1 The nuclear amplitude 



The nuclear breakup amplitude is obtained from eq.(|ll|) by expanding to first 
order in V e f / and separating the term which is zero order in V e f / . We assume 



that the interaction time is so short that the w-dependence in eq.(ll) can be 
neglected. This is the sudden approximation or 'frozen halo approximation'. 
The integral can be evaluated by changing the time variable to t' = t + z/v. 
The nuclear breakup amplitude reduces to 

CT'(k,d) = yd 3 re- lkr (exp(-z Xnt (b))-l)^ m (r) (13) 
d 2 r±e - tk± .r ± (exp (_ iXnt (b)) _ 1} {r ^ kz) 

where the eikonal phase 
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Xnt (b) = \ [ V nt (rx+d, t) dt. (14) 

The vector b = + d is the impact parameter of the neutron relative to the 

target. The breakup amplitude ( |i~3| ) is equivalent to the one used in |^(| [l3| 
except that in those references the dependence of the neutron-target interaction 
on the combination z + vt is taken into account before setting u> = 0. The effect 
is to replace the z-component k z of the final momentum by k\ where 

h = k z + - = k z + (15) 



3.2 The Coulomb amplitude 

The calculation of Coulomb breakup is different from nuclear breakup because: 
i) The Coulomb force has a long range, ii) The Coulomb force does not act 
directly on the neutron but it affects it only indirectly by causing the recoil of 
the charged core. The recoil effect is included in the effective interaction V e //. 
The Coulomb contribution can be obtained by putting the Coulomb potential 
in the definition (||) and is 



V«//(r,R(i)) = + A ZpZ T e 2 (?^5Wj (16) 

= +AW (17) 

[cP + {vt) 2 ) 

Equation (jl^) is the effective Coulomb potential on the neutron in the dipole 
approximation. The effective force on the neutron (|5|) is in the opposite direction 
to R(t). 

The fact that the Coulomb interaction has a long range means that the de- 
pendence of the breakup amplitude on the excitation energy e^ — Eq = hcu cannot 
be neglected. In the present section we give an expression for the breakup am- 
plitude which is correct to first order in the Coulomb interaction and which 
also contains the effects of the neutron-target interaction as well as the in- 
dependence. The Coulomb contribution to the breakup amplitude gf^ ul (k,d) 
is obtained from jTT| ) by writing V2 = V n t + V e f / and expanding to first order 
in V e f f . The first order contribution to the integrand in (|llj) is proportional to 

e -ix»«(r,t) (y eff ( r , t) _ %// ( r>t) y nt ( r , t) ) 

= V eff (r, t) - (1 - e -*»«(r,t)) V eff (r, t) - ie-^^Xeff (r,t) V nt (r, i)(18) 

where 

Xnt (r,t) = i / V nt (r, t') dt', xef f (r,t) = i / V eff (r, t 1 ) dt 1 (19) 



1 f 00 1 f 

; (r,t) = - y (r, < Xe// (r,t) = -j 
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The qualitative effects of the neutron-target interaction in the Coulomb 
breakup amplitude for any value of u> can be understood by making a simple 
approximation to the nuclear part of the eikonal phase in eq. (Eq) ; 



X»t(r,t)=Xnt(b)0(t), (20) 

where 9 (t) = 0, t > 0; 8 (t) = 1, t < 0. The physical assumption is that the 
neutron target interaction is important only for a short time around t = corre- 
sponding to the point of closest approach between the projectile and the target. 
The Coulomb effective interaction is more slowly varying with time and has a 
contribution for t > when the nuclear phase eq. ( pfj| ) is zero and a contribution 
for t < when the nuclear phase is important. With this approximation the 
Coulomb amplitude can be written as a sum of two terms 

g?™ 1 (k,d) = (M) + 9& (k,d) (21) 



The first term is obtained from the first term in eq.(|18[). It depends only 
on the Coulomb interaction and yields the standard perturbation result for 
Coulomb breakup in the dipole approximation 

gZ l (k,d) = -if d 3 re-* * Xeff (r, W ) 0, m (r) (22) 



where 



1 f°° 

Xe//(r,u;) = - J dte«*V eff (T,t) (23) 

The second part of ( |2~]f ) has contributions from the second and the third 
terms in eq. ( |T^ ) . In the third term we treat again the neutron-target interaction 
in the sudden approximation. Then the sum of these two contributions is 

(k,d) = J d 3 re- 4k ' r (l - exp {-i Xn t (b)))B (x, z, w) 4> Un (r) (24) 



nC 
9lm 



where 



i 

B (x, z, w) = — 
n 



dte tut V eff (T,t)+ / dtV eff (r,t) 



(25) 



Explicit expressions for g^ ert (k,d) and B (x, z, w) are given in the appendix. 
The amplitude (|24| ) has an interesting structure. It is exactly like the nuclear 
breakup amplitude but with an effective neutron wave function 

cj>efj (r, vj) = B (x, z, w) 0o (r) (26) 

If 4>o (r) is an s-state then <fi e ff (r, w) is a p-state. 

The results presented in the next section have been obtained with the pre- 
vious formulae, where the sudden approximation or frozen halo approximation 
has been made only for the neutron-target interaction. The same approxima- 
tion can also be made for the Coulomb amplitude by putting the Coulomb 
adiabaticity parameter w — cod/v equal to zero in the above expressions. When 
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huj = e k - e fa 1 MeV, E lab = 41 A.MeV and d « R s = 11.5 fm then to pa 0.10. 
This number is relatively small so the sudden approximation might also be rea- 
sonable for the Coulomb amplitude in some cases. Thus the expression for the 
complete amplitude (|ll|) in the sudden approximation becomes 

g lm (k,d) = J <f > r± e- ik± - r± (exp (-i{ X nt (b) + Xeff (b))) - 1) 4> im {r ± ,k z ) 

(27) 

The Coulomb part to first order is 

gf£* (k,d) = -i / d 2 r ± e- lk — Xe // (b) exp (-t Xnt (b)) 4> lm (r Xj *,) (28) 
which is consistent with eqs. (IDII) and (||. 
3.3 Cross sections 

This section contains the formulas needed for calculating neutron breakup cross 
sections. The final state for the 3-body system in the breakup reaction is spec- 
ified by the momentum k conjugate to coordinate r and the momentum K 
conjugate to R. They are related to final momenta of the core, neutron and 
target by 

k c = -k + £ 2 K, k n = k + /3xK, k T = -K. 

In the eikonal approximation for the core-target scattering the full 3-body 
breakup amplitude A (K, k) is given by 

A (K, k) = J d 2 de- lK ^ d S ct (d) g (k, d) (29) 

where S c t (d) is the profile function describing the scattering of the halo core 
by the target. It is given in terms of the core-target potential by the eikonal 
formula 

S ct (d) = exp (~ J dtVct (d, vt)j . (30) 
The differential cross-section integrated over Kj_ is 

Jt = ^/ d2K±|A(K ' k)|2 = ^3/ d2d l^( d )| 2 |3(k,d)| 2 . 

If j3\ is small then k„ is almost the same as k. If k f = (k x , k y) kt) is the neutron 
final momentum in the target reference frame, while k = (k x , k y , k z ) is in the 
projectile reference frame and k z = kf z — m n v/h, = k 2 , + ky, then the above 
cross section is equivalent to the expression for the cross section for emitting 
neutrons with energy Sf into the solid angle dVL. It is obtained by integrating 
the breakup probability distribution over the impact parameter d of the center 
of mass of the projectile relative to the target p8[ 
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where P e i = \S c t (d) | 2 is the probability that the projectile core-target system 
remains in the ground state during the reaction. Several simple parameteriza- 
tions are possible for P e i- The profile function S c t (d) can be calculated from 
the core-target optical potential with the eikonal formula (^). It can also be 
approximated by the strong absorption model 

P e i = 1 for d > R s 

P e i = for d < R s (32) 

where R s is the strong absorption radius for the core-target collision. In this 
way the calculations are analytical up to the cross section formula ^6|. The 
breakup probability is obtained from 



d 2 P{d) _ 1 m n kf 1 \ auucI | A Coui\2 /qo\ 

le~7dn~8^ +A lm |. (33) 

J m 



The initial neutron bound state wave function in the projectile has angular 
momentum (I, m) and the amplitudes Af^f 1 and Af^ ul are related to gi m (k,d) 
by a phase factor 

A i m + Am = e fa (k,d) . (34) 

The relation between the neutron-core relative energy and the neutron en- 
ergy in the laboratory is 

e k = Sf + -jm n v 2 - kf z vti. (35) 

Eq. (p4|) makes a connection with the notation in earlier works j68| . The 
phase factor has its origin in a different choice of the coordinate system. 

In cq.(j33|) the momentum of the breakup neutron has polar angles (0,4>) 
where the azimuthal angle <$> refers to the theoretical reaction plane defined by 
the incoming beam momentum kj n and the impact parameter d. This plane is 
not an observable but when K±d >> 1 then it is almost the same as the plane 
defined by the incoming direction and the direction of K = k c + k„. This is 
because the main contribution to the integral ( p9|) for A (K, k) comes from the 
range of angles where is almost parallel to or almost antiparallel to d. 



4 Results 

In this section we discuss some sample results for the breakup of the halo state 
in 11 Be obtained by using the formalism of the previous sections, and in order 
to asses the accuracy of the formalism we will compare to existing experimen- 
tal data. Initial state parameters are the same as in p6fl . Namely the initial 
neutron binding energy for the 2s state is Eq = — 0.5MeV and the spectroscopic 
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Table 1: Energy dependent optical model parameters for ?i+ 9 Be. a#=0.387fm, 
r,=1.368fm, a/=0.3fm at all energies. R R = (1.447A^ /3 - 0.005(e/ - 20))/m. 





V R 


AW S 


W v 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


20-40 


38.5-0. 1456/ 


1.666+0.365e/ 


0.3756/ - 7.5 


40-120 




16.226-0.1(e/ -40) 


7.5-0.02(6/ - 40) 


120-180 




8.226-0.07(6/- 120) 


5.9 



factor C 2 S — 0.77 [Q. The neutron-target optical potential, denned as in 
[jB7| , for the Be target is given in Table 1 and it is the same as |27j for the 
other targets. We will show separately the cross sections due only to the nu- 
clear amplitude, to the Coulomb amplitude and their coherent and incoherent 
sums. The Coulomb amplitude eq.(21) is a sum of two parts: the perturba- 
tion amplitude eq. (E2) which depends only on the Coulomb interaction and the 



Coulomb- nuclear term eq. ( 24 ) which contains the distorting effect of the nuclear 
interaction on Coulomb breakup. Both neutron energy spectra and angular dis- 
tributions will be discussed. 

The main purpose of the calculations presented here is first to clarify some 
differences in the mechanisms of Coulomb and nuclear breakup, then to asses 
the importance of the Coulomb-nuclear coupling term and of the interference 
effects. Their influence on the final neutron energy distribution is discussed 
in section 4.1 and on the energy-integrated angular distribution in section 4.2. 
The first results for the angular distribution in section 4.2 are integrated over 
the azimuthal angle of the outgoing neutron direction. This corresponds to an 
experimental situation where the final direction of the halo core is not measured. 

In the second part of the discussion in section 4.2 the angular distributions 
are integrated over a limited range of azimuthal angles measured with respect 
to the plane of the classical trajectory. As explained at the end of the first part 
of section 3 this plane is almost the same as the plane defined by the incoming 
direction and the direction of K = k c + k„ provided |Kj_| is large enough. In 
any case a measurement of such an angular distribution requires a measurement 
of both k„ and k c . 



4.1 Energy distributions 

In Figs, (la) and (lb), we have represented the nuclear and Coulomb breakup 
spectra for different values of the angle 9 as indicated, in the case of the Au 
target. In Fig. (lb) the Coulomb cross section at very forward angles ( 6 — 0) is 
zero for e/ w m n v 2 /2 and it has then two asymmetric peaks. The double hump 
structure in the Coulomb breakup cross-section is a known feature of dipolc 
Coulomb breakup from DWBA (eg. calculations by Okamura |54|). However it 
was stated in |j4j that the origin of the double peak was not clear in the DWBA 
formalism. Its origin is instead very clear in our semiclassical treatment. 
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Figure 1: Neutron energy distributions following nuclear (a) and Coulomb (b) 
breakup of n Be at 4lA.MeV at fixed 6 ( 197 Au target). Solide line 6 = degree, 
dashed line 6 = 1.5 degrees, doted line 6 = 3 degrees, dashed-dot line 6 = 6 
degrees. 



In fact from the approximation used in Eq.(^) one sees that the core target 
interaction gives rise to an effective interaction on the neutron Eq. (|^) . Then as 
Coulomb breakup is a reaction which transfers some momentum to the neutron, 
there cannot be breakup when the incident neutron keeps both its incident ve- 
locity and direction. This kinematical effect can be easily understood by looking 
at Eqs.(|22|) and (^). The initial wave function momentum distribution has a 
maximum at k z = for s and p states. In such conditions kf z = kf-cosd = ^f^- 
and the neutron going at zero degrees should keep all the momentum available 
from the relative motion. In fact the Coulomb operator contains a derivative 
with respect to k z which then gives a zero probability for such a process. On 
the other hand if the neutron is slightly deflected it can keep his incident en- 
ergy with a maximum probability. The peaks at the other angles correspond to 
kinematical situations where the neutrons have all the available incident energy 
per nucleon at the distance of closest approach d corresponding to the velocity 
v = (2(Ecm — VcB{d))/ [i) 1 / 2 where fj, is the projectile-target reduced mass. 
Such an energy is slightly less than the incident energy per particle because 
of the Coulomb barrier. This is always the case for nuclear breakup as shown 
in Fig. (la). The kinematics of Coulomb breakup can be easily understood 
by looking at Fig. (2) where we show the relevant momentum vectors in the 
laboratory reference frames, as indicated. The incident momentum per particle 
is mv and p', p" are the two neutron momenta in the projectile reference frame 
corresponding to the same given neutron laboratory angle 6. |p'| and |p"| are 
equal to the averaged momentum transferred to the neutron by the Coulomb 
effective field. Along the laboratory final momentum axis we sketch the cor- 
responding energy distribution as already given in Fig. (lb). It is clear that 
increasing the angle 6 the two peaks tend to join each other, merging finally to 
just one peak at a critical angle 6 C = y/ 2mVc b /Ap / fiv beyond which Coulomb 
breakup is kincmatically suppressed. Note that 6 C is the angle at which the 
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Coulomb breakup angular distribution is a maximum, as shown later in Fig. (5). 




Figure 2: Kinematics of the breakup reaction. The solid curve represents the 
Coulomb energy distribution for a given small 8. 

On the other hand, in the nuclear breakup, Fig. (la), the neutron average 
momentum in the projectile reference frame is k z w 0, corresponding to the 
maximum of the initial momentum distribution and of the breakup probability. 
Then the energy distribution has just one peak at all angles, close to the incident 
energy (and momentum) per particle. 

Double peaked longitudinal distributions are typical of reactions in which 
the mechanism of breakup involves an effective repulsive force between the frag- 
ments. For example in heavy-ion fragmentation reactions the so called Coulomb 
rings are observed |7l| . This happens because the reaction goes through a kind 
of inelastic excitation in the projectile which then decays in flight. Therefore 
it is suggested that exclusive measurements of longitudinal distributions with 
the neutron in coincidence with the core, for fixed, small, angles 9 could help 
distinguishing different reaction mechanisms. This could be very useful in the 
case of two-neutron halo breakup in which the second neutron decays in flight 
[ p9[ |30| , |3l| from a resonant state and therefore with a mechanism different from 
the nuclear breakup which is responsible for the decay of the first neutron. 

We show in Figs. (3) the neutron final energy spectra for Coulomb and for 
nuclear breakup for the reactions 9 Be( 11 Be, 10 Be + n) 9 Be, 48 Ti( 11 i?e, 10 Be + 
n) i8 Ti, 197 Aii( n i?e, 10 Be + n) 197 Au at £ mc =4lA.MeV §5| in the target ref- 
erence frame. We have used the same low energy cutoff (27MeV) as in the 
experiment. We sec that the nuclear breakup (dot-dashed line) is dominant for 
the 9 Be target, for the 48 Ti target Coulomb (long dashed line) and nuclear give 
similar contributions, and for 197 Au, Coulomb breakup becomes dominant. The 
Coulomb-nuclear coupling term (dotted line) increases with the target mass and 
for 197 Au is of the same order as the pure nuclear term. From these results, 
we would expect big interferences between Coulomb and nuclear processes for 
48 Ti and 197 Au. In fact the solid line which represents the coherent sum of the 
three processes is rather different from the dotted line with crosses, which is the 
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Figure 3: Neutron energy distributions following breakup of 11 Be at 41A.MeV 
for several targets ( 9 Be, 48 Ti, 197 Au). Nuclear breakup distributions are repre- 
sented by the dot-dashed line, Coulomb by the long dashed line, the Coulomb- 
nuclear coupling term by the dotted line, their coherent sum by the solid line. 
The dotted line with crosses is the nuclear plus the perturbative Coulomb inco- 
herent sum. 



incoherent sum of nuclear and perturbative Coulomb breakup. 

Finally in Fig. (4a) we show the neutron energy spectrum in the projectile 
reference frame together with the data from Q for Coulomb breakup on 208 Pb 
at 72A.MeV. The same spectrum is also given in the target reference frame in 
Fig. (4b). Notation is the same as in Fig.(l). In this case as in the case of the Au 
target discussed before the contribution of the Coulomb-nuclear coupling term 
is very close to that of the pure nuclear term. Our results agree well with the 
experimental spectrum. We used the known spectroscopic factor C 2 S = 0.77 
[ [72| for the s-initial state in 11 Be. The integrated Coulomb cross section is 
1.70b to be compared with the experimental value of 1.8± 0.4b. This is to show 
that our treatment of the Coulomb breakup by perturbation theory is accurate 
within the experimental error bars for the description of the shape and absolute 
magnitude of Coulomb breakup. It has been suggested p|, [l^, |ll| |2(], |2l| that 
a treatment beyond perturbation theory could be more accurate in same cases. 
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Proton breakup, in particular at lower energies than those discussed here, is 
one of such cases. On the other hand methods that have studied consistently 
neutron breakup p], [l0| in first order or higher order, within the same model, 
have found negligible differences. Non perturbative treatments are also expected 
to be more accurate in the case of smaller neutron separation energies and lower 
beam energies j2(j] than those discussed here. The spectrum of Fig. (4) has been 
calculated and discussed in a number of theoretical works [|l0| [lij [l5| which all 
use different non-perturbative methods but whose results are all consistent with 
each other and with our results, within the experimental error bars. 
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Figure 4: Neutron energy distributions following breakup of n Be on 208 Pb at 
72A.MeV in the projectile frame (a), and in the target frame (b). The notation 
is the same as in Fig.(l). Experimental points are from Nakamura [4]. 



4.2 Interference effects in the neutron angular distribu- 
tions 

Fig. (5) contains the data from the exclusive experiment of Ref. [^5) in which the 
neutron angular distribution following the breakup from 11 Be was measured in 
coincidence with the 10 Be ejectile. Three different targets were used to check on 
the relative importance of the nuclear and Coulomb breakup mechanisms. In the 
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Figure 5: Neutron angular distributions following breakup of 11 Be at 41A.MeV 
for several targets ( 9 Be, 48 Ti, 197 Au). In the bottom figures the nuclear distri- 
bution is represented by the solid line, Coulomb by the long dashed line, and 
the Coulomb-nuclear term by the dotted line. In the top figures the coher- 
ent sum of the nuclear, Coulomb and Coulomb-nuclear terms by the solid line, 
the nuclear plus the perturbative Coulomb incoherent sum by the dashed line. 
Experimental points are from Anne p5[ . 



bottom part of the figure the solid line is the nuclear breakup cross section, the 
long dashed line is the Coulomb breakup while the dotted line is the Coulomb- 
nuclear coupling term. In the top part of the figure we show the experimental 
data together with the coherent sum of the three contributions (solid line). The 
dashed line is the nuclear plus the perturbative Coulomb incoherent sum. It is 
important to notice that the large angle scattering is due mainly to the nuclear 
breakup and it is sensitive to the neutron-target optical potential used. About 
two thirds of the total amount of nuclear breakup is due to large angle scattering. 
Therefore new measurements of neutron angular distributions, including large 
angle scattering data would greatly help in settling the question of the relative 
importance of nuclear and Coulomb breakup for heavy targets. On the other 
hand spectra of the type shown in Fig. (4a) correspond to small angle scattering 
only (typically 6 < 15deg) where the nuclear is very small and therefore fitting 
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Figure 6: Neutron angular distribution as in Fig. 5, but 0< <f> <30 degrees for 
several targets ( 9 Be, 48 Ti, 197 Au). The results in the top figures are integrated 
over the core target impact parameter d. The bottom figures are results at a 
fixed impact parameter d = R s . The solid line represents the coherent sum of 
nuclear, Coulomb and Coulomb-nuclear terms. Their incoherent sum by the 
dashed line. 

The theoretical results for the neutron angular distributions were summed 
over the neutron energy spectrum and averaged over its out-of-plane angle <fr 
because this is what is contained in the experimental data. The peak values of 
the cross section show clearly an increase with the target due to the Coulomb 
breakup. From the comparison of the coherent and incoherent sums it is clear 
that coupling and interference effects arc more important at the angles where 
nuclear and Coulomb terms are closer in magnitude, which happens at small 9 
for the 9 Be target and around 20-30deg for the other two targets. 

We have seen that the effects due to the interference are partially washed out 
by the average over <p. On the other hand we have identified the reason of the 
lost of coherence at large in the integration of the probabilities over the core- 



target impact parameter d as contained in Eq.(31). To clarify these points we 
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show in Fig. (6) similar calculations for the same targets as before but this time 
we have chosen to restrict the (j> integration to the range -j- 30deg. In the top 
part the integral over the core-target impact parameter has been performed. In 
the bottom part of the figure the calculations have been performed integrating 
one fermi out of the strong absorption radius R s . The results of these figures 
are however still summed over the neutron energy spectrum. 




Figure 7: Cosinus of the Coulomb-nuclear phase as a function of the impact 
parameter for several 8 angle and et — 4lMeV (fixed k z , k y = 0, k± = k x ). 
The solid line stands for 9 = 6 degrees, the dashed line for 9 = 12 degrees and 
the dot-dashed line for 9 — 18 degrees. The target is 197 Au. 

The solid lines are the results for the cross section obtained by summing 
coherently the nuclear and Coulomb and Coulomb-nuclear amplitudes while 
the dashed line is again the incoherent sum of the nuclear plus the perturbative 
Coulomb. From the bottom figure it is clear that the interferences change sign 
as a function of 9. 

Eqs.(|l6|) and ( p2| ) show that the important contributions to the integral 
for the perturbative Coulomb amplitude come from the projectile region while 
integral ( |l~3| ) for the nuclear amplitude is weighted towards the target. As a 
consequence the Coulomb-nuclear phase difference has a dependence on the 
projectile-target impact parameter. The cosine of the relative Coulomb-nuclear 
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Table 2: Integrated (9 = -j- 60deg) breakup cross sections in barn. Incident 
energies in A.MeV. 



Target 


E- 

^inc 


Coul-Nucl 


Nuclear 


Coulomb 


Incoh. sum 


Cohe. sum 


% 


y Be 


30 


0.6 10~ a 


0.200 


0.011 


0.212 


0.222 


5 




41 


0.5 10~ 3 


0.191 


0.009 


0.200 


0.209 


5 




120 


0.4 10~ 4 


0.044 


0.004 


0.048 


0.050 


5 


4« xi 


30 


0.2 10" 1 


0.260 


0.292 


0.573 


0.621 


9 




41 


0.2 10" 1 


0.252 


0.229 


0.496 


0.531 


8 




120 


0.6 10" 2 


0.234 


0.098 


0.338 


0.374 


10 


197 Au 


30 


0.327 


0.372 


2.98 


3.68 


3.50 


-5 




41 


0.245 


0.354 


2.36 


2.96 


2.85 


-4 




120 


0.151 


0.321 


1.05 


1.52 


1.53 


1 




72 


0.148 


0.346 


1.70 


2.19 


2.17 


-1 



phase as a function of d is shown in Fig. (7) for a set of neutron emission angles. 
The neutron is emitted in the plane defined by the impact parameter and the 
incident beam direction. The cosine of the relative phase has a regular oscillatory 
behavior and the period of the oscillations decreases with increasing 9. The 
oscillations reduce the magnitude of the Coulomb-nuclear interference. They 
also influence its sign as indicated by the integrated cross sections in Table 2. 

Finally we give in Table 2 the absolute cross sections for each individual term: 
the Coulomb-nuclear amplitude, the eikonal nuclear and Coulomb perturbative. 
Also their coherent and incoherent sums are shown, for the targets studied in 
this work. We have integrated in the experimental range 9 = -j- 60deg. In 
the calculations of this work we have used the strong absorption hypothesis 
Eq.(p2[) with the prescription R s = 1.4(Ap 3 + A^ 3 )fm. Using the smooth 
absorption hypothesis |6(| would reduce the absolute nuclear cross section but 
leave the energy and momentum distributions unchanged ||. For example the 
total coherent cross sections on 9 Be at the three given energies, with a smooth 
cut off are 0.221b, 0.200b and 0.042 b respectively. The same would happen 
using the new formalism proposed in [ fL3| . The values without interference are 
very similar to those already obtained by some of us in [ gTJ . We remind the 
reader that when the neutron is detected in coincidence with the core only the 
diffraction-type of nuclear breakup contributes. Cross section values are given 
at three incident energies E inc = 30, 41 and 120A.MeV and at 72A.MeV for the 
208 Pb target. In the last column we give the percentage of interference effect. 
The effect of the interference on the angle integrated energy distribution is of 
the order of 4-j- 10% therefore small as already seen by other authors || [ll], p2| . 
The interference is constructive for the 9 Be and 48 Ti targets but destructive 
for the 197 Au and 208 Pb targets. For the heavy targets it appears also that the 
interferences tend to vanish at high energy. Then from the analysis of the energy 
and angular distributions we conclude that the effect of the new Coulomb- 
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nuclear coupling term and of the interferences is to give a slight increase with 
respect to the simple sum of the nuclear plus the Coulomb perturbative breakup. 
However from the analysis of the results for the total cross sections it appears 
that the total interference of the three terms can be destructive with respect to 
their sum for the heavy target. This is manly due to the imaginary part of the 
neutron target optical potential. The situation might be different in the case of 
the proton breakup studied in 17 24|, |TJ |2^] because the effective charge in the 
Coulomb potential Eq.(|36|) will be different from ours obtained for a neutron. 

At this stage it is difficult to draw a general conclusion on the effects of 
the interference and of the Coulomb-nuclear coupling as they clearly depend on 
the target and on the incident energy and on the angular range covered by the 
experiment. 

The results discussed in this paper are restricted to the 2s initial state in 
11 Be. The Coulomb breakup from the d-componcnt of the ground state gives 
a contribution two orders of magnitude smaller than the s-component and thus 
it has been neglected here. However the formulae given in the appendix hold 
for any initial angular momentum. Calculations for cases like the weakly bound 
carbon isotopes, where k = 2 orbitals are important, will be presented elsewhere. 



5 Conclusions 

In this work we have derived an eikonal formalism leading to Eq.(p"l|) to treat 
consistently the nuclear and Coulomb breakup to all orders in the interactions. 
This formula could be studied numerically. Here we have studied only the effect 
of the first order coupling and interference between the two processes. By taking 
first order terms in the Coulomb potential and the nuclear potential to all orders 
we have obtained a scheme in which the breakup amplitude is a sum of three 
terms. Two of them reduce to the well known forms of the nuclear breakup 
in the eikonal model and of the Coulomb breakup in first order perturbation 
theory. The third term, Eq. ( gjj) which is derived here for the first time, depends 
on both the Coulomb and nuclear potentials and can be viewed as one of the 
lowest " higher order coupling terms" . Within such a model numerical results 
show that interferences between nuclear and Coulomb breakup are present in 
halo nucleus reactions and in the data already existing. They are responsible 
for the peak at small angles in the angular distribution on the 9 Be target and 
for the behavior at large 8 in the case of heavier targets. They have not been 
seen so far more clearly in the data because of the averaging over the neutron 
out-of-the-reaction plane angle <f) and the integration over the core-target impact 
parameter. However they would show up in neutron double differential angular 
distribution from halo breakup by making a selection on a small range of <j) 
angles and detecting the core at fixed angles, provided it follows a Coulomb 
trajectory. In this case the reaction plane could be identified and a one by 
one correspondence between the core detection angle and the trajectory impact 
parameter could be applied. 

The Coulomb-nuclear interference effects are small for the cases studied in 
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the this paper. Interference terms do not have a definite sign in the angular 
distribution but they are overall constructive in the energy spectra and absolute 
cross sections at medium energies and for light to medium mass targets. For 
heavy targets the results presented here show an overall destructive interference 
which depends on the imaginary part of the neutron target optical potential 
and which tends to vanish increasing the beam energy. The effect of the new 
Coulomb-nuclear coupling term is rather small and negligible for light targets, 
while for heavy targets is of the same order as the pure nuclear term. This is 
evident in the energy distributions and in the small angle angular distributions. 
The comparison with the data and with other works is reassuring in so far as 
the neglect of higher order terms in the Coulomb amplitude is concerned. They 
scale as powers of the effective charge, that is as powers of C = PiZpZre 2 /hv 
and f3i is small (/3i = 0.1) for the 11 Be projectile. Higher order terms will be 
very small for light targets but could become more important for heavy targets. 
They have not been studied here but they can be evaluated within the adiabatic 
approximation Eq. ( p7| ) . 

Thus it appears that the so-called Coulomb experiments would need a more 
accurate theoretical analysis than performed so far. Also we think that it would 
be very important to measure more neutron angular distributions, like the ones 
[ p5[ discussed in this paper. They are a beautiful tool to identify the reaction 
mechanism and to test the theoretical methods. Therefore we suggest that in 
more exclusive kind of experiments it should be possible to see better the effects 
discussed in this work. It will be necessary to have more intense radioactive 
beams giving also a better statistics and to use neutron position sensitive de- 
tectors over a wide angular range. New facilities of the type being planned in 
Europe (EURISOL) and USA (RIA) are therefore most welcome. 
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A Details of the Coulomb and Coulomb-nuclear 
amplitudes. 

The Coulomb contribution to the effective interaction is 



We need to evaluate several integrals. The first is the Fourier transform 



V eff (r,t) = +p 1 Z P Z T e 2 



xd + zvt 



(36) 



(d? + (vi) 2 ) 



3/2 




(37) 



By making a change of variables s — vt/d the first integral becomes 



21 



(l + s 2 ) 3/2 
C 

= 2— {xvjKi (zu) + izzuKo (tu)) (38) 

where 

v nv 

and the constant C = PiZpZtz 2 /hv is a dimensionless interaction strength. 
The functions^ an d K\ are modified Bcssel functions and can be defined by 
the integrals 



1 f cos (rut) dt 

2 ./_«, (l+<2)3/2 



„ . , I f°° tsm(zat)dt 

Substituting eq.pS) into eq.(p2|) yields the perturbation amplitude 



of r 

gP ert (k,d) = -i— / d 3 re- lk r (xwK x (w) + izzuK Q (zu)) $ lm (r) 

= 2- [wK l {m)-^-+imK (s {m)^- ) tm (k). (42) 

The constant C = fiiZpZ^e 2 /Tiv is a dimensionless interaction strength and 
is the adiabaticity parameter. The functions -FTo and ifi are modified Bessel 
functions and <fii m (k) is the Fourier transform of cf>i m (r). First order pertur- 
bation theory in V e ff should be a good approximation if CR p /d < 1/2 where 
Bp is the radius of the halo nucleus. Numerical values for the modified Bessel 
functions show that the adiabatic or 'frozen halo' approximation should be sat- 
isfactory if vj < 0.2. Eq.(|4^) is the standard dipole expression for Coulomb 
breakup. 

The function 



i 

B (x, z,vj) = — 
n 



(43) 



dte wt V eff {r,t)+ / dtV eff {r,t) 

involves two integrals. The first is 

1 f°° C 

- j dtV eff (r,t) = -(x + z) (44) 
The other is a one-sided Fourier transform 

-j_ dte™ t V eff (r,t) = -J dte-«*V eff (r,-t) (45) 
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It can be expressed in terms of zuKq (zu) , zuK\ (zu) and two new integrals 



f> , ^ f°° sin (zut)dt 

wKi(a) - /„ irhk (46) 

- , , f°° tcos(zut)dt , A . 

" A ° < " ) - L ~o^?f T (47) 

Putting everything together 

C 

B (x, z, zu) = i — ix(wK\ (zu) + 1 — izuK\ (zu)) — z (zuKq (zu) — 1 — izuKo (zu))) 
a 

(48) 

B Initial states with general angular momenta. 
B.l The Coulomb breakup amplitude. 

If the initial state wave function is approximated by its asymptotic form which 
is an Hankel function j73j 

^m(r) = -i l C ilo hf ) {i"ior)Yi m {e,<}>), 7o r » 1, (49) 



where Ci is the asymptotic normalization constant |7Cj ] and 70 = \] — 2meo/h, 
then the general form of the initial state momentum distribution is given by the 
Fourier transform of Eq. ( flg| ) 

~4> lm {^k z ) = ^ x -^—^— Yum {k) (50) 

where k = (k x ,k y ,k z ) is a real vector. 

In nuclear breakup reactions the use of the asymptotic part of the initial 
bound state wave function is justified because of the short range nature of the 
nuclear interaction which gives breakup form factors in Eq. ( |l3| ) localized in the 
overlap region between the two interacting nuclei. The amplitude for Coulomb 
breakup is expressed in terms of the Fourier transform of the bound state wave 
function in Eq . (f42| ) . This Fourier transform is well approximated by the Fourier 
transform Eq.(p0) of the corresponding Hankel function provided that 70-R << 1 
and kR << 1 where R is the nuclear radius. These conditions are often satisfied 
in the Coulomb breakup of a halo nucleus. We have checked this for 11 Be using 
wave functions calculated in a square-well potential and have shown that for 
|eq| < lMeV and up to I = 2 the contribution of the internal part of the wave 
function to the full Fourier transform can be neglected. 

We want to evaluate 2n < J2 m \ Af^ ul | 2 > ip where <> v stands for the av- 
erage over the <p angle of ky. Some definitions for the neutron momenta are: 
k/ = (k x . k y . kfj) is the neutron final momentum in the target reference frame, 
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while k = (k x , k v ,k z ) is in the projectile reference frame and k z = kj 

;„2 _ l2 , ju2 
_L x ~~ ^ y ' 

We define 



A 



Coul 



F /fn (k) 



where ifi and are modified Bessel functions and 
Flrn(k) = fi(k)Y lm (k) 

k l 



fi(k) = 



7^(7o 2 



with 70 = y/—2meo/h. Then 

B^Y^FU^FUk) 



21 + 1 
4tt 



f(k')f(k)Pi(cose) 



i n v/u, 

(51) 

(52) 
(53) 

(54) 



with cos9 = k.k'/fc'fe. 

We calculate the derivatives like d 2 BJ dk' x dk x and finally put k' = k. Note 
that <9/dfc'Pi(cos60| k=k , = 0. 

Using the above results in Eq. (|5l|) and averaging over <p 



2n<Y,\A?r l \ 2 >v 



21 + 1 



v 2(k\ + klf 

+ K 2 (co)[f' 2 k 2 k 2 z /2 + P[(l)f 2 kl]) 
After simplifications, the final result is 



(Kf(u,) [f 2 k 2 k 2 j2 + P{(l)f(kl + kj2)} 



(55) 



2n<J2\A^ l \ 2 > v = 



21 + 1 



K\ [Ui + -k\yi] + Kl [Hi + klVi] ) (56) 



with 



U t {k) 



Vi(k) = 



1(1 + 1) 



u2i-2 



2 j 2l (j 2 + k 2 ) 2 
k 21 - 4 [1(1-1) 



llKll + k 2 ) 4 



(57) 



(58) 



B.2 The Coulomb cross section. 

We give in the following the analytical formulae for the cross section due to 
Coulomb breakup only obtained using 

Eqs.d^HIf). In the integral Eq.© 
over the impact parameter d we change variable of integration from d to w 
given in Eq. ( |39|) . Furthermore to compare to experimental angular distributions 
we average over the angle <j>. The expressions for the double differential cross 
section, in the case of an initial s and p state are 



d a Coui mkf > 2 



n(kl+k 2 z+1 2 ) 4 



[A(K 2 - Kl) + BK Ki\ (59) 
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where 
and 

B = 4^L (61) 

CO 

in the case of an s-initial state, while for a p-initial state we get 

A = i (ki + (kl + 7o 2 ) 2 + 2{k ^ +1 l) [kl + ll ~ klf + Akl{2kl + 7 ° 2) ) 

(62) 

B = (| (*1 + (kl + 7o 2 ) 2 + 55^g) <* + 7o 2 - «)")) (63) 
ui> — uj{d — R s ) is the argument of the modified Bessel functions in Eq.(S). 
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